Method for Receiver Autonomous Integrity Monitoring and Fault 
Detection and Elimination 

FIELD OF THE INVENTION 

[0001] The present invention relates generally to Fault Detection and Elimination 
(FDE) in a discrete-time controlled process, and particularly to methods for Receiver 
Autonomous Integrity Monitoring (RAIM) in global positioning systems (GPS). 

BACKGROUND OF THE INVENTION 

[0002] GPS uses satellites in space to locate objects on earth. With GPS, signals from 
the satellites arrive at a GPS receiver and are used to determine the position of the 
GPS receiver. Currently, two types of GPS measurements corresponding to each 
correlator channel with a locked GPS satellite signal are available for civilian GPS 
receivers. The two types of GPS measurements are pseudorange, and integrated 
carrier phase for two carrier signals, LI and L2, with frequencies of 1.5754 GHz and 
1.2276 GHz, or wavelengths of 0.1903 m and 0.2442 m, respectively. The 
pseudorange measurement (or code measurement) is a basic GPS observable that all 
types of GPS receivers can make. It utilizes the C/A or P codes modulated onto the 
carrier signals. The measurement records the apparent time taken for the relevant code 
to travel from the satellite to the receiver, i.e., the time the signal arrives at the 
receiver according to the receiver clock minus the time the signal left the satellite 
according to the satellite clock. The carrier phase measurement is obtained by 
integrating a reconstructed carrier of the signal as it arrives at the receiver. Thus, the 
carrier phase measurement is also a measure of a transit time difference as determined 
by the time the signal left the satellite according to the satellite clock and the time it 
arrives at the receiver according to the receiver clock. However, because an initial 
number of whole cycles in transit between the satellite and the receiver when the 
receiver starts tracking the carrier phase of the signal is usually not known, the transit 
time difference may be in error by multiple carrier cycles, i.e., there is a whole-cycle 
ambiguity in the carrier phase measurement. 

[0003] With the GPS measurements available, the range or distance between a GPS 
receiver and each of a multitude of satellites is calculated by multiplying a signal's 
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travel time by the speed of light. These ranges are usually referred to as pseudoranges 
(false ranges) because the receiver clock generally has a significant time error which 
causes a common bias in the measured range. This common bias from receiver clock 
error is solved for along with the position coordinates of the receiver as part of the 
normal navigation computation. Various other factors can also lead to errors or noise 
in the calculated range, including ephemeris error, satellite clock timing error, 
atmospheric effects, receiver noise and multipath error. With standalone GPS 
navigation, where a user with a GPS receiver obtains code and/or carrier-phase ranges 
with respect to a plurality of satellites in view, without consulting with any reference 
station, the user is very limited in ways to reduce the errors or noises in the ranges. 

[0004] To eliminate or reduce these errors, differential operations are typically used 
in GPS applications. Differential GPS (DGPS) operations typically involve a base 
reference GPS receiver, a user GPS receiver, and a communication mechanism 
between the user and reference receivers. The reference receiver is placed at a known 
location and the known position is used to generate corrections associated with some 
or all of the above error factors. The corrections are supplied to the user receiver and 
the user receiver then uses the corrections to appropriately correct its computed 
position. The corrections can be in the form of corrections to the reference receiver 
position determined at the reference site or in the form of corrections to the specific 
GPS satellite clock and/or orbit. Corrections to the reference receiver position are not 
as flexible as GPS satellite clock or orbit corrections because, for optimum accuracy, 
they require that the same satellites be observed by the user receiver and the reference 
receiver. 

[0005J The fundamental concept of Differential GPS (DGPS) is to take advantage of 
the spatial and temporal correlations of the errors inherent in the GPS measurements 
to cancel the noise factors in the pseudorange and/or carrier phase measurements 
resulting from these error factors. However, while the GPS satellite clock timing 
error, which appears as a bias on the pseudorange or carrier phase measurement, is 
perfectly correlated between the reference receiver and the user receiver, most of the 
other error factors are either not correlated or the correlation diminishes in wide-area 
applications, i.e., when the distance between the reference and user receivers becomes 
large. 
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[0006] To overcome the inaccuracy of the DGPS system in wide-area applications, 
various wide area DGPS (WADGPS) techniques have been developed. The 
WADGPS includes a network of multiple reference stations in communication with a 
computational center or hub. Error corrections are computed at the hub based upon 
the known locations of the reference stations and the measurements taken by them. 
The computed error corrections are then transmitted to users via a communication link 
such as satellite, phone, or radio. By using multiple reference stations, WADGPS 
provides more accurate estimates of the error corrections. 

[0007] Thus, a user with a GPS receiver may use different modes of navigation, i.e., 
standalone GPS, DGPS, WADGPS, carrier-phase DGPS, etc. Whichever of the 
navigation modes is used, there is always the possibility that the range with respect to 
a satellite are computed based on a faulty measurement, such as a measurement with 
respect to a failed satellite. When this range is used in determining the position of the 
user, an erroneous or wrong position would result. Thus, a faulty measurement can 
cause serious degradation to the reliability and integrity of the GPS system. Therefore, 
various integrity monitoring techniques have been developed for fault detection and 
elimination (FDE) in GPS systems. Receiver autonomous integrity monitoring 
(RAIM) is the name coined by the FAA for methods of integrity monitoring in GPS 
using redundant GPS satellite measurements. 

[0008] The literature on RAIM and FDE procedures is extensive. Most of the 
procedures in the literature, however, are related to aviation use and attempt to bound 
the probable error in a position domain. As a result, they generally involve very 
extensive computations. One of the earliest papers describing a RAIM procedure is a 
paper by Brown and McBurney, "Self-Contained GPS Integrity Check Using 
Maximum Solution Separation," Navigation, Vol. 35, No. 1, pp 41-53. In this paper, 
the authors suggest: (1) obtaining GPS measurements with respect to n satellites in 
view; (2) for each of the n satellites, solving for the user position based on 
measurements with respect to the other (n-1) satellites; (3) computing all possible 
distances between the solutions in the horizontal plane and determining a maximum 
distance among the possible distances; and (4) using the maximum distance as a test 
statistic and declaring a failure when the maximum distance exceeds a threshold. 
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Clearly, this technique is very computationally intensive and does not isolate a 
particular measurement or satellite as being faulty. 

[0009] Another early paper is by Parkinson and Axelrad, "Autonomous GPS Integrity 
Monitoring Using the Pseudorange Residual/' Navigation, Vol 35, No. 2, pp 255-271. 
In this paper, the authors suggest an excellent test statistic based upon pseudorange 
measurement residuals, but when it comes to using the test statistic to isolate a failed 
satellite, they use a scheme similar to that used by Brown and McBurney, i.e., for 
each of a plurality of satellites, they compute a test statistic while leaving out the 
measurement with respect to the satellite. Again, this procedure presents an excessive 
computational burden. 

SUMMARY 

[0010] A method for detecting and identifying a faulty measurement among a 
plurality of GPS measurements, obtained by a GPS receiver with respect to a plurality 
of satellites in view of the GPS receiver, determines whether the plurality of GPS 
measurements include a faulty measurement. In response to a determination that the 
plurality of GPS measurements include a faulty measurement, the method identifies a 
satellite contributing the faulty measurement by computing a correlation value 
associated with each of the plurality of satellites, and selecting a satellite associated 
with a highest correlation value as the satellite contributing the faulty measurement. 
In one embodiment, to insure that the correct satellite is identified, the satellite 
associated with the highest correlation value is selected when the highest correlation 
value exceeds a predetermined threshold value and the predetermined threshold value 
is sufficiently larger than a second highest correlation value. In an alternative 
embodiment, the satellite associated with the highest correlation value is selected 
when the difference between the highest correlation value and a second highest 
correlation value exceeds a predetermined threshold. 

[0011] In some embodiments, whether the GPS measurements include a faulty 
measurement is determined by computing test statistic using post- fix residuals 
corresponding to the plurality of GPS measurements, and comparing test statistic with 
a threshold residual value, which is chosen based on a navigation mode used by the 
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GPS receiver. If the test statistic exceeds the threshold residual value, a faulty 
measurement is detected. 

[0012] In some embodiments, the correlation value associated with a satellite is the 
absolute value of a correlation coefficient associated with the satellite. The correlation 
coefficient is computed based on a residual sensitivity matrix corresponding to the 
plurality of satellites and a residual vector including the post-fix residuals 
corresponding to the plurality of GPS measurements. 

[0013] In some embodiments, the size of the error in the faulty GPS measurement is 
determined based on a residual sensitivity matrix corresponding to the plurality of 
satellites and the root mean square residual. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] Figure 1 is a block diagram of a computer system that can be used to carry out 
a method for detecting and identifying a faulty GPS measurement or a satellite 
contributing to the faulty GPS measurement. 

[0015] Figure 2 is a flowchart illustrating the method for detecting and identifying a 
faulty GPS measurement or a satellite contributing to the faulty GPS measurement. 

[0016] Figure 3 is a flowchart illustrating a method for determining whether a 
plurality of GPS measurements include a faulty GPS measurement. 

[0017] Figure 4 is a flowchart illustrating a method for identifying a faulty GPS 
measurement among a plurality of GPS measurements. 

[0018] Figure 5 A is a flowchart illustrating a method for identifying a satellite with a 
highest correlation value as the satellite contributing to the faulty GPS measurement. 

[0019] Figure 5B is a flowchart illustrating another method for identifying a satellite 
with a highest correlation value as the satellite contributing to the faulty GPS 
measurement. 

[0020] Figure 6 is a flowchart illustrating a method for verifying that the satellite 
contributing to the faulty measurement has been identified correctly. 
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DESCRIPTION OF EMBODIMENTS 



[0021] FIG. 1 illustrates a computer system 100 that can be used to carry out the 
method for detecting and identifying a faulty GPS measurement among a plurality of 
GPS measurements. Each of the plurality of GPS measurements is taken by a GPS 
receiver 122 based on signals from one of a plurality of satellites 1 10-1, 1 10-2, . . ., 
1 10-n, where n is the number of satellites in view of the GPS receiver 122. The 
plurality of satellites, or any one or more of them, are sometimes referred to hereafter 
in this document as satellite(s) 1 10. In some embodiments, the GPS receiver 122 and 
the computer system 100 are integrated into a single device, within a single housing, 
such as a portable, handheld or even wearable position tracking device, or a vehicle- 
mounted or otherwise mobile positioning and/or navigation system. In other 
embodiments, the GPS receiver 122 and the computer system 100 are not integrated 
into a single device. 

[0022] As shown in FIG. 1, the computer system 100 includes a central processing 
unit (CPU) 126, memory 128, an input port 134 and an output port 136, and 
(optionally) a user interface 138, coupled to each other by one or more 
communication buses 129. The memory 128 may include high-speed random access 
memory and may include nonvolatile mass storage, such as one or more magnetic 
disk storage devices. The memory 128 preferably stores an operating system 131, a 
database 133, and GPS application procedures 135. The GPS application procedures 
may include procedures 137 for implementing the method for detecting and 
identifying the faulty GPS measurement, as described in more detail below. The 
operating system 131 and application programs and procedures 135 and 137 stored in 
memory 128 are for execution by the CPU 126 of the computer system 124. The 
memory 128 preferably also stores data structures used during execution of the GPS 
application procedures 135 and 137, including GPS pseudorange and/or carrier-phase 
measurements 139, as well as other data structures discussed in this document. 

[0023] The input port 134 is for receiving data from the GPS receiver 122, and output 
port 136 is used for outputting data and/or calculation results. Data and calculation 
results may also be shown on a display device of the user interface 138. 
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[0024] FIG. 2 illustrates a method 200 for detecting and identifying a faulty GPS 
measurement among a plurality of GPS measurements obtained by the GPS receiver 
122 with respect to the plurality of satellites 1 10. As shown in FIG. 2, method 200 
includes step 210 for determining whether the plurality of GPS measurements include 
a faulty measurement. In response to a determination in step 210 that the plurality of 
GPS measurements include a faulty measurement, method 200 further includes step 
220 in which the faulty measurement is isolated or identified among the plurality of 
GPS measurements, or the satellite contributing to the faulty measurement is isolated 
or identified among the plurality of satellites. With the satellite contributing the fault 
measurement identified, method 200 may include an optional step 230 in which a size 
of an error in the faulty measurement is determined, and an optional step 240 for 
verifying that the correct identification has been made. 

[0025] FIG. 3 illustrates an embodiment of a method 300 for determining whether the 
GPS measurements include a faulty measurement in step 210. As shown in FIG. 3, 
method 300 includes step 310 in which post-fix residuals corresponding to the 
plurality of GPS measurements are computed and step 320 in which a test statistics is 
formed using the post-fix residuals and is compared with a threshold to determine 
whether the plurality of GPS measurements include a faulty measurement. A residual 
of a GPS measurement represents a disagreement between the GPS measurement and 
a prediction or expected value of the GPS measurement. Before the position and clock 
bias of the GPS receiver are adjusted, the residuals are often referred to as pre-fix 
residuals or measurement innovations. A measurement innovation can be computed 
based on the difference between a GPS measurement and a theoretical prediction of 
the GPS measurement. Alternatively, a measurement innovation corresponding to a 
GPS measurement can be computed as the difference between the GPS measurement 
and an expected value of the GPS measurement computed from an initial estimated 
state of the GPS receiver, as discussed in the following. 

[0026] Whatever the navigation mode, navigating with GPS involves a discrete-time 
controlled process that is governed by a linear stochastic difference equation: 

Hx = z + n (1) 
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where a: is a state vector of the discrete-time controlled process and, in the case of 
GPS, it includes corrections to the position and clock bias of the GPS receiver, H is a 
measurement sensitivity matrix including direction cosines of the state vector or unit 
vectors from the GPS receiver to each of the n satellites, z is a measurement 
innovation vector including measurement innovations corresponding to the plurality 
of GPS measurements, and n is a measurement noise vector. In the case of carrier- 
phase measurements, the state vector may include an unknown ambiguity factor. 

The post-fix residuals are usually obtained in a two-step process. First, a least-squares 
solution for x is made, i.e. 

x = (H T Hy l H T z (2) 



or, 



x = (H T R- l Hy l H T R~ l z 



(3) 



for a weighted least squares solution. R in Equation (4) is a measurement covariance 
matrix and 



R = 



0 



0 



(4) 



where cr s , i = 1 , 2, . . n, represents a standard deviation of GPS measurement noises 

with respect to the z th satellite. An example of the methods for calculating cr. can be 

found in "Precision, Cross Correlation, and Time Correlation of GPS Phase and 
Code Observations,^ by Peter Bona, GPS Solutions, Vol. 4, No. 2, Fall 2000, p. 3-13, 
or in "Tightly Integrated Attitude Determination Methods for Low-Cost Inertial 
Navigation : Two-Antenna GPS and GPS/Magnetometer" by Yang, Y., Ph.D. 
Dissertation, Dept. of Electrical Engineering, University of California, Riverside, CA 
June 2001, both of which are hereby incorporated by reference. 

[0027] The correction to the state as calculated in Equation (3) or (4) is used to 
transform the measurement innovations (pre-fix residuals) into a set of post-fix 
residuals in step 310 of method 300, according to the following equations 
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(5) 



where A is a residual vector including as its elements post-fix residuals corresponding 
to the plurality of GPS measurements, and S is a residual sensitivity matrix: 

S = {\-H(H T HY X H T ) (5a) 

or, 

S = (1 - H(H r R- l HY l H T R' 1 ) (5b) 

[0028] S is called the residual sensitivity matrix because it is a matrix whose elements 
are the residuals corresponding to unit changes in the measurement innovations. This 
can be explained through the following discussions. Equation (1) can be expanded to 
include a set of state vectors corresponding to a set of arbitrary measurement 
innovation vectors: 

HX = Z (6) 

where ^includes a set of state vectors jci,jc2, ...,x n corresponding to a set of arbitrary 
measurement innovation vectors zi,Z2,... f Zn.« Now, if we let Z be an identity matrix, we 
get the state vectors corresponding to a set of measurement innovation vectors each 
representing a single unit change (e.g., one meter or any arbitrary unit) in the 
measurement innovation of a different respective one of the plurality of satellites 
while the measurement innovations of the other satellites are kept unchanged, as 
shown in the following equation: 

HX = I (7) 
The least-squares solution for X is then: 

X = (H T Hy l H T (8) 
For a weighted least-squares solution this becomes: 

X = (H T R~ l H)~ l H T R~ l (9) 
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[0029] Multiplying equation (8) by //gives a prediction about the measurement 
innovations of each satellite based on the least-squares solution: 



HX = H(H T H)~ l H T 



(10) 



Subtracting this prediction from the input value of the innovations, i.e., the identity 
matrix, we obtain a matrix including the residuals of unit changes in the measurement 
innovations, which is the residual sensitivity matrix S in Equation (5a) or (5b). As 
shown in Equation (5), the S matrix can be used to directly map the pre- fix residuals 
(measurement innovations) into the post-fix residuals. 

[0030] Thus, each column (or row) of the S matrix includes residuals corresponding 
to a unit change in the measurement innovation of one of the plurality of satellites 
while the measurement innovations of the other ones of the plurality of satellites are 
kept unchanged. The S matrix has a number of interesting properties: it is symmetric; 
it is idempotent, i.e., S Z =S 2 =S 3 =... ; the sum of the elements in any row or column 
equals zero; and the length of any row or column is equal to the square root of an 
associated diagonal element. Since the state vector, x y has four elements for most 
navigation modes, the rank of S is n-4 t where n is the number of satellites. 

[0031] With the S matrix, the post-fix residuals in the residual vector A can be 
computed according to Equation (5) in step 310. The post-fix residuals can be used to 
compute a root mean square residual 5, which is the norm value (or length) of a 
residual vector A divided by the number of measurements, n, i.e., 



[0032] After computing the post-fix residuals in step 310, method 300 further 
includes step 320 in which the post-fix residuals or the root mean square residual are 
used to form a test statistic, which is then compared with a fault threshold to 
determine whether the plurality of GPS measurements include a faulty measurement. 
The test statistic can be the root mean square residual. Or, it can be the length of the 
post-fix residual vector scaled by appropriate normalization for the number of 
satellites, such as the square root of (n-4), as in the following equation 




(11) 
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(»-4) 



(12) 



Alternatively, when one is not concerned with very small errors, the fault threshold 
can be set large enough so that the scaling is relatively insignificant. Thus, the test 
statistic a can simply be the length of the post-fix residual vector, i.e., 



[0033] When the test statistic is larger than the fault threshold, it is determined that 
the plurality of GPS measurements include a faulty measurement. Since the level of 
measurement noise (or the position accuracy) of the GPS receiver 122 often depends 
on the navigation mode in which the GPS receiver is operating, in some embodiments 
the fault threshold value is selected to correspond to the navigation mode used by the 
GPS receiver 122. 

[0034] FIG. 4 illustrates an embodiment of a method 400 for identifying the faulty 
measurement or a satellite contributing the faulty measurement in step 220 of method 
200. As shown in FIG. 4, method 400 includes step 410 in which a correlation value 
associated with each of the plurality of measurements or with each of the plurality of 
satellites is computed, and step 420 in which a satellite associated with a highest 
correlation value among the plurality of satellites 1 10 is identified as the satellite 
contributing the faulty measurement, or the GPS measurement with respect to the 
satellite associated with the highest correlation value is identified as the faulty 
measurement. 

[0035] In some embodiments, the correlation value associated with each one of the 
plurality of satellites is computed based on a correlation coefficient associated with 
the satellite or with the one of the plurality of GPS measurements with respect to the 
satellite. The correlation coefficient associated with a satellite, e.g., the j th satellite, 
where j = 1 , 2, . . ., n, represents a correlation between the residuals of the plurality of 
GPS measurements and the residuals corresponding to a unit change in the 
measurement innovations with respect to the j th satellite, while the measurement 
innovations of the other satellites are kept unchanged. Thus, the correlation 



(12a). 
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coefficient associated with the j satellite can be computed using the S matrix and the 
residual vector A according to the following equation: 

' mJ w 

where p j is the correlation coefficient associated with the j th satellite, A. is the i th 
element of the residual vector A , sf is the element in the i th row and j th column of the 
S matrix, |a| is the norm value of the residual vector A, s J represents the j th column 
of the S matrix, and the summation is over all elements of A or of the j th column of 
the S matrix. Note that the length of sf 9 \s j \ , is equal to the square root of the diagonal 

element, sj . 

[0036] In one embodiment, the correlation value associated with a satellite is equal to 
the absolute value of the correlation coefficient associated with the satellite. With the 
correlation values thus computed, the satellite associated with a highest correlation 
value among the correlation values associated with the plurality of satellites can 
generally be identified in step 420 as the satellite contributing to the faulty 
measurement. However, an additional check may be required to ensure that the 
correct satellite is identified. 

[0037] When only five satellites are available, the correlation coefficient associated 
with every satellite should be about one (or minus one). This is because the degree of 
freedom is only one (one more satellite than the number of unknowns in the state 
vector). With only five measurements corresponding to the five satellites, each 
column in the S matrix is perfectly correlated with the residual vector, except that the 
length of the columns are different. In other words, the large root mean square 
residual, which indicated a fault in step 210 of method 200, could be caused by a 
faulty measurement from any one of the satellites. Moreover, even when the number 
of satellites is greater than 5, a similar phenomenon can occur with a degenerate 
geometry. Such a degenerate geometry is indicated if more than one satellite are 
associated with correlation coefficients close to one or minus one. 
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[0038J Thus, to insure that the satellite contributing to the faulty measurement be 
identified correctly, method 500A or 500B, as shown in FIGS. 5A or 5B, respectively, 
can be used in step 420. Method 500A includes step 51 OA in which the highest 
correlation value is identified among the correlation values associated with the 
plurality of satellites. Method 500A further includes step 520A for determining that 
the highest correlation value exceed a predetermined threshold value and step 5 3 OA 
for determining that the predetermined threshold value is sufficiently larger than a 
second highest correlation value among the correlation values associated with the 
plurality of satellites. Alternatively, method 500B is used, which includes step 51 OB 
in which the highest and the second highest correlation values are identified, step 
520B in which the difference between the highest and second highest correlation 
values are computed, and step 530B for determining that the difference exceeds a 
predetermined minimum difference value. 

[0039] With the satellite contributing the faulty measurement identified, method 200 
may further include step 230 in which the size of any error in the faulty measurement 
(or the size of the fault) is estimated. Estimating the size of the fault is sometimes 
useful, especially when method 200 is used to identify a faulty measurement in a 
Real-Time Kinematic (RTK) navigation mode. Since RTK navigations usually 
involve carrier phase measurements and thus resolutions of whole-cycle ambiguities, 
a measurement fault can be the result of a cycle slip in a tracking loop or an improper 
determination of the whole-cycle ambiguity. When this is the case, the size of the 
fault will be a multiple of the carrier wavelength. Given that the satellite, such as 
satellite 1 10-k, where k is 1,2, . . ., or n, with the highest correlation value is 
identified, a best estimate of the size of the error, e k , is given by the root mean square 
residual divided by the length of the column of the S matrix associated with the 
satellite 100-j. So, is the root mean square residual corresponding to a unit change 
in the measurement innovation of satellite-k while the measurement innovations for 
the rest of the plurality of satellites are kept unchanged. Since the length of the 
column of the S matrix is equal to the square root of the corresponding diagonal 
element of the S matrix, we have: 



(14) 
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Note that with poor geometry the diagonal element can be small. However, the 
smaller the diagonal element, the larger the measurement error from the 
corresponding satellite must be before it can cause the root mean square residual to 
exceed the detection threshold. 

[0040] With the size of the error in the faulty measurement estimated, method 200 
may further include step 240 for verifying that there is only one faulty measurement 
among the plurality of GPS measurements and that the identification of the satellite 
contributing the faulty measurement is correctly made. FIG. 6 illustrates a method 600 
for performing the verification in step 240. As shown in FIG. 6, method 600 includes 
step 610 in which the residual vector A is adjusted to account for the error according 
to the following equation: 

A' = A + SE (15) 

where E is a vector with zero elements corresponding to satellites whose 
measurements are without error and an element having the value of for the satellite 
100-k that contributed the faulty measurement. Method 600 further includes step 620 
in which the test statistic is recalculated according to Equation (1 1), (12), or (12a) 
after the adjusted residual vector A' is used to replace the residual vector A in the 
equation. If the fault identification has been made correctly, the recalculated test 
statistic should now pass the threshold test for the test statistic, indicating that the 
fault has been removed. Thus, method 600 further includes step 630 in which the 
recalculated test statistic is compared with the fault threshold to verify that there is 
only one faulty measurement among the plurality of GPS measurements and the 
identification of the satellite contributing the faulty measurement is correctly made. 

[0041] The above described method for identifying a faulty GPS measurement is 
computationally efficient because it works entirely in the measurement domain. When 
there are more than five satellites being tracked, the method can be used to 
dramatically improve the reliability of positioning and navigation using GPS, with 
only a small cost in additional computational complexity. 
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